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We examine the role of spin twists in the formation of domain walls, often called stripes, 
, by focusing on the spin textures found in the cluster spin glass phases of La2-xSrxCu04 and 

Yi-xCaxBa2Cu306. To this end, we derive an analytic expression for the spin distortions pro- 
duced by a frustrating bond, both near the core region of the bond and in the far field, and then 
derive an expression for interaction energies between such bonds. We critique our analytical theory 
by comparison to numerical solutions of this problem and find excellent agreement. By looking at 
collections of small numbers of such bonds localized in some region of a lattice, we demonstrate the 
stability of small "clusters" of spins, each cluster having its own orientation of its antiferromagnetic 
order parameter. Then, we display a domain wall corresponding to spin twists between clusters of 
^ , locally ordered spins showing how spin twists can serve as a mechanism for stripe formation. Since 

the charges are localized in this model, we emphasize that these domain walls are produced in a 
situation for which no kinetic energy is present in the problem. 
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I. INTRODUCTION: 



''O ■ The now frequent experimental observations of spin and/or charge modulations in the cuprate superconductors 
Ch and related doped transition metal oxides [Q was predicted by the "frustrated phase separation" phenomenology of 
^ , Emery and Kivelson. Their theoretical considerations involved the assertion that a doped Mott insulator phase 
separates as a consequence of the competition between the kinetic energy of mobile holes and the magnetic energies of 
an antiferromagnetic (AFM) phase with long-range correlations. While it seems unlikely that this scenario is correct 
in the strong correlation limit when one adds the Coulomb interaction into this problem Emery and Kivelson 
argued that macroscopic phase separation became "frustrated" , and the resulting anomalous normal state possessed 
low-energy fluctuations corresponding to stripes, or domain walls Q. To connote that these entities correspond to 

■ metallic stripes, it is now common to refer to these structures as "rivers of charge" 
Recent neutron scattering studies |^ of the single-layer La2-z;Sr2,Cu04 (LSCO) system have revealed that (at least 

in experimental results to date) an elastic magnetic response associated with static stripe-like correlations are only 
found in (i) low x systems (x < 0.06) at low temperatures such that the transport is that of a doped semiconductor [0, 

■ and (ii) in the famous a; « 1/8 system. One could interpret both the strongly disordered, low temperature results, and 
the a; ~ 1/8 data, as evidence that pinning effects are necessary, either from disorder or commensurability interactions, 
to produce static stripes. Such arguments are consistent with the successful approach of Tranquada and coworkers in 
producing static stripes that could then be observed in scattering experiments in a variety of systems [Q . 

. The low X, low temperature region of the LSCO phase diagram in which the static magnetic stripe correlations are 

■ found corresponds to the spin-glass phase of LSCO. Early magnetic resonance work on this system |^ suggested that 
Q , it is appropriate to think of this phase as a cluster spin glass, so named because small clusters of spins achieve their 
O ■ own short-range AFM correlations, but the cluster-cluster ordering is spin-glass like. A numerical simulation of this 
^ ' phase, coupled with new crystals and new susceptibility data, lended support to this characterization |^ . 

In the latter paper |^ , one conundrum associated with the mechanism behind the formation of the cluster spin glass 
phase was pointed out, and goes as follows: The frustrated phase separation phenomenology claims that support for 
such physics is found in the existence of the cluster spin glass phase Q]. However, detailed analysis of the transport 
in this region of the LSCO phase diagram concludes that the transport is similar to that of a doped semiconductor 
[0. Thus, at least in the low-temperature cluster spin-glass phase the competition between kinetic and magnetic 
energies does not exist in the form proposed originally by Emery and Kivelson, and the question can then be asked, 
does the absence of the holes' kinetic energy from extended states not eviscerate the frustrated phase separation 
phenomenology as a viable mechanism associated with the formation of the cluster spin glass phase? Put another 
way, the numerical simulations of Rcf. that found evidence for the character of the spin-glass phase being like that 
of a cluster spin glass produced this spin texture with zero kinetic energy, and thus is the kinetic energy a necessary 
ingredient in the formation of stripe phases pO[|? 
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This question becomes more important in view of recent experiments of Julien, et al. These NMR/NQR 

resuhs demonstrated the existence of a so-cahed charge glass at higher temperatures, foUowed by the appearance of 
the superconducting phase at lower temperatures, followed by the cluster spin glass phase at the lowest temperatures. 
The frustrated phase separation phenomenology predicts this sequence of charge glass/cluster spin glass phases, and 
thus unlike the above arguments, suggests that the cluster spin glass is stabilized by the pinning of the charge stripes 
by defects, followed by the subsequent freezing of the spin degrees of freedom within clusters defined by the pinned 
stripes of the charge glass phase. Unfortunately, again, the transport of this system at low temperatures is insulating 
(for, say, T < 75 K), so in the spin glass phase there are no rivers of charge which could "carve out" the domain walls 
of the cluster spin glass phase! 

In this manuscript we present analytical results that formalizes the claim that via quenched disorder from the Sr 
impurities one can produce the topology of pinned stripes without any kinetic energy. To this end, we derive the 
spin distortion pattern produced by such quenched disorder (which frustrates the background AFM order), and then 
demonstrate that this theory successfully predicts the stability of localized clusters of AFM correlated spins produced in 
a situation with zero kinetic energy. In this case, the origin of the stripes associated with the domain walls comes from 
the spin twists between the clusters, the clusters themselves having been produced by spin twists of the spin texture 
as the background spins attempt to accommodate the frustrating magnetic interactions produced by the quenched 
disorder. A well known extrapolation to higher temperatures [ [l2| then implies that the qualitatively identical spin 
twists generated by mobile carriers must be part of the mechanism associated with the formation of stripes. 

We wish to make clear that our paper does not claim to be the first to propose that magnetic interactions in 
general, and spin twists in particular, are important in the formation of stripes. Firstly, the work of Salem and one 
of us [|3[ investigated the problem of frustrating FM bonds whose locations could be chosen such that the ground- 
state energy was minimized. It was found that when quantum fluctuations were included, if the magnitude of the 
frustrating interaction was smaller than that of the background majority spins, periodic stripes of frustrating bonds 
were the ground state configuration. (So, in these ground states, again there is no kinetic energy hut stripe phases 
are indeed encountered.) Secondly, when the frustrating bonds cannot chose their (static) positions but are fixed by 
the Sr impurity ions, the numerical simulations (mentioned above) of Ref. [|| suggested the presence of domain walls 
between the clusters of the cluster spin-glass phase. More recently, work of Stojkovic and coworkers examined 
a version of the mobile hole problem by implementing a purely magnetic model that included the long-ranged spin 
twists produced by mobile holes (as well as the frustrating Coulombic energy between the carriers) and found many of 
the magnetic structures encountered in ||l3|| , including stripe phases. Lastly, White and Scalapino, who find evidence 
for stripe structures in the t-J model for mobile holes note that the charge and spin distributions of striped 
structures attempt to accommodate the frustration (read: spin twists) on the magnetic background produced by the 
mobile holes. 

Looking at the totality of the evidence in this and the above-mentioned papers we believe that one can make a 
strong case that there is a similar mechanism at work in the formation of stripes in all of these situations, and that 
this mechanism is spin twists. 

Our paper is organized as follows. In the next section we present a detailed theoretical analysis of the effects of 
quenched disorder on the spin texture in systems such as weakly doped LSCO, producing a reliable analytical theory 
of the spin distortions both near the frustrating bond and in the far field. We use this distortion field to produce 
an accurate interaction functional between pairs of such bonds, which we then use to demonstrate the stability of 
such clusters in the cluster spins glass phase. In particular, this leads to a clear identification of the local AFM order 
parameter of each cluster. Finally, we show the resulting spin texture between two such clusters, and demonstrate 
how (local) stripe configurations can be stabilized in the cluster spin-glass phase. 



We consider the familiar model of magnetism in the CuO planes of the high cuprates in which the copper ions 
and oxygen holes are treated as a lattice of spins governed by a Heisenberg Hamiltonian 



where {ij) denotes a summation over near neighbour pairs of spins and the Jij are the exchange interaction integrals. 
For an undoped lattice these spins are the Cu spins and the exchange integrals are equal and negative; the ground 
state corresponds to an AFM ordered state on a square lattice. In what follows we simplify our considerations by 



II. CORE SOLUTION AND ENERGIES OF THE SINGLE BOND PROBLEM: 



A. Hamiltonian and definitions 




(1) 
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using classical spins, implying that only the transverse (viz., moment reorientation) and not the longitudinal {viz., 
moment magnitude) spin-spin interactions are included |12| . 

The simplest model of the effects of doping in weakly doped cuprates at low temperatures was proposed by Emery, 
and corresponds to localizing the holes on oxygen sites and replacing the AFM Cu-Cu superexchange for this occupied 
bond with an effective FM exchange. The phase diagram of the multiply doped version of this model was produced 
by Aharony, et al . Although detailed transport analysis of this part of the LaSrCuO phase diagram has shown 
that a slightly different model (l^J^ of the localized dopants is required for a direct comparison to experiments, the 
FM bond model (which we shall refer to as the frustrating bond model) is more amenable to analytical study, for 
reasons that we shall elaborate on below, and shall be used throughout this paper. 

Thus, we consider the Hamiltonian of Eq. (]^) wherein the i,j label sites of a square lattice (that is, only the Cu 
spins and the effective interactions between them are considered) and the exchange interaction integral between two 
adjacent sites i and j has the form 

J _ ( XJ with probability x/2 , . 

•-^'^ ^ I - J with probability l~x/2 

where J and A are positive constants and A represents the relative strength of the ferromagnetic and antiferromagnetic 
bonds jl^. The doping level x could be, say, either the Sr doping level in La2_xSrxCu04 or the Ca doping level in 
Yi_xCaxBa2Cu306 (noting that Neidermayer et al. jl^] has shown that the phase diagrams in these two systems are 
identical). If we now choose a coordinate system such that linear combinations of x- and y-directed unit vectors span 
the lattice, then the Hamiltonian can be written explicitly as 

H = — — Ji^i-^-aSi ■ Si^a (3) 

i a 

where i is summed over all lattice sites and a ranges over ±x, ±y. The equilibrium condition corresponds to that of 
zero torque from the local effective field at each lattice site: 

^ J,,,+aS, X S,+a -0. (4) 

d 

As is well known, the complication of treating a bipartite lattice (labelling the two sublattices as A and B sites) 
can be avoided by transforming the physical problem of FM bonds in a predominantly AFM background into the 
mathematically equivalent problem of AFM bonds in a FM background. These two pictures can be converted one to 
the other under the simple transformation (AFM — > FM) given below: 

Jij ^ (5) 

+S, for i e A 
-S, for i e B ■ 

Lastly, we note that the objects under consideration are classical spins, and we set their length to be one, scaling J 
to be J . Further, the ground states that we discuss in this paper all correspond to situations in which the spins lie 
in some plane, and thus from now on we restrict our formalism to describe planar spins. We denote the bulk direction 
of the spins by Sqo , and at any lattice site i there exists a spin characterized by the angle between the spin and 
the 2:-axis: 

Si = X cos ijji + y sinV'i ^ • Sj = cos {ipi — ipj) . (6) 

We choose ipi = + V'oo where ipoo is taken to be the average angle of the spins over the bulk of the material. That 
is Soo — xcostAoo + ysinipoo so that the angle 0^, defined according to cos (f>i — Si ■ Sao, represents the deviation of 
the spin at site i from the bulk direction. The collection {(pi} of the spin distortions at each lattice site constitutes 
the spin texture on the lattice. 

The numbering scheme for the lattice sites near the frustrating bond (what from now on we call the bond sites) is 
shown in Fig. |^. 
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FIG. 1. The lattice site numbering scheme. The frustrating bond sits between the and 1 sites. 



B. Spin deviations of a frustrating bond 

In a FM lattice doped with a single AFM bond, the ground state solution to the spin texture is no longer obvious, 
and such a situation is called frustrated. We wish to produce an accurate analytical solution to this problem in both 
the far-field region and close to the frustrating bond. This will allow us to accurately track the energies of both the 
single and many bond problems. 

Consider a purely FM lattice frustrated by the introduction of an x-directed AFM bond between the (0, 0) and 
(1,0) lattice sites. We shall denote the spin distortions at these sites by (po and 0i, respectively. The Hamiltonian is 
then written as 

El 
cos {(f>i — (/'j) + A J cos {(f>o — (f>i) (7) 

(y) 

where the prime indicates the omission of the AFM bond from the summation. It is no longer clear that the trivial 
solution (f>i = represents the ground state since there may be another solution with tt/2 < \(j>o — (j>i\ < 3tt/2 that the 
takes the system to a lower energy state. 

As is well known, considerable information can be gained from an examination of the dynamical properties of 
a linearized system of equations of motion about a proposed equilibrium structure. Here, we briefly outline this 
formalism, since it will be important to our later work. For our purposes, the dynamical behaviour of a lattice of n 
independent spins is modeled sufficiently by 

-^'i = - ^ sin (0i - (/)i+a) j = l,2,...,n. (8) 

a 

Close to the ordered equilibrium state = 0, this behaviour is governed by the linearized system x ~ — Afx where 
M = Df{0) is the derivative matrix of / evaluated at the origin and x is related to the spin texture according to the 
row vector x^ = (/)2, 03, . . . , ^„]. 

Solving for the normal modes, and taking note of negative eigenvalues, the instability of the system to a non FM 
ordered ground state can be identified. This technique, when applied for larger and larger systems, reproduces the 
known result pO[ | that an instability is first reached at Ac = 1 and that there is only one stable spin texture for all A 
exceeding Ac. 

Unlike such instability analysis, or the Fourier-based approach of Vanninemus, et al. pO| , here we wish to develop 
a continuum theory capable of describing an infinite lattice including the spins in the immediate neighbourhood of the 
frustrating bond. To this end, we proceed as follows. 

Let (/) he a, function of a continuous variable r which ranges over the entire xy-plane such that i— > (j){vi). Then, 
provided that is a smooth, slowly varying function of position, we may approximate the equilibrium condition for 
the undoped system (to lowest order [£1[) by 

V2(/. = 0. (9) 
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Now consider a FM lattice frustrated by the introduction of a single a;-directed AFM bond. The equilibrium 
spin distortions away from the core of the frustration are governed by Laplace's equation. Choose the origin of the 
coordinate system centred on the bond, and solve Laplace's equation, in polar coordinates, by separation of variables. 
The imposition of the appropriate solution symmetries yields 

oo 

(p{r,e)= r^'^An.cosme (10) 

m— 1,3,5,... 

where we have adopted the convention that summations are over odd indices only. 

It is clear that for sufficiently large r, the lowest order term dominates. That is to say, far from the bond the 
distortions are dipolar: 

0(r) = (11) 



where p = Aix or p = Aiy. This agrees with the well known results in the literature (see, e.g., Refs. p(]| , pi| ). 
Unfortunately this result is inadequate for our purposes, since we also require the spin distortions near the bond. As 
we show below, a previous attempt |2|] fails, and thus we present new arguments leading to a valid solution close to 
and far away from the frustrating bond. Other work has been unable to solve analytically this part of the frustrating 
bond problem (although it is clear that they are aware of the issues that we have finally solved). 

We have written down a general solution to the static spin texture on the infinite lattice due to a single AFM bond, 
and that solution consists of a linear combination of an infinite number of possible solution modes. However, since we 
have shown that the single bond system has only one stable solution, we expect that any prepared state will decay 



into the state of lowest energy given by the m = 1 solution in Eq. (10). That is every for m ^ 1. Nonetheless, 

we run into the difficulty that the field equation to which Eq. ( p^ ) is a solution is not strictly valid at the bond sites. 
Consequently, we cannot expect that these solutions will perform well near the bond itself. Indeed, we find that the 
purely dipolar solution with 1/r fall-off fits numerical solutions extremely well (~ 0.5%) up to three or four lattice 
sites away from the bond, but that in the core of the frustration the deviation becomes quite large. At the bond sites 
themselves, the error is ~ 25%. 

We may ask, of course, why such a description does not suffice if, for the most part, we are interested in the spin 
distortions away from the bond. Surely we can tolerate a small error at a handful of lattice sites? The answer is that 
we cannot. Since the spin distortions are most severe in the immediate vicinity of the bond, the spins in the core of 
the frustration are a large contributor to the magnitude of the total energy stored in the distortions. Thus a proper 
calculation of the energy in the system requires that we model the core correctly. Equally important is that, for a 
given solution mode, the local equilibrium condition at the bond site determines the overall magnitude of the spin 
distortions. That is, it fixes the magnitude of the dipole moment associated with the distortion field. Kovalev and 
Bogdan, who suggested a continuum approach for the core region of this problem , fall into precisely this trap and 
hence obtain the wrong magnitude for the long range behaviour for the spin distortions (see below). 

It remains to be answered how we might treat the spin distortions around the bond. Ideally, we would like to 
treat the AFM character of the interaction in the core as a small perturbation on the field equations, but this is not 
possible since the continuum formalism is badly behaved at the origin under the symmetries we have imposed. A 
second possibility would be to derive a separate discrete solution valid in the core and to match it smoothly onto the 
exterior continuum solution. However, we are inclined to avoid such a patch-work approach. Not only is it somewhat 
inelegant, but it also defeats the purpose of introducing the continuum formalism, namely, to do away with discrete 
calculations altogether. 

Instead, we make use of the fact that the sets {cos{m9)/r™} and { sin (m6')/r'"} are complete (in the sense that any 
static spin texture satisfying the symmetries [ p2| specified by a single x- or y-directed bond can be expanded in one 
of these bases). Thus, to solve for the spin distortions everywhere on the lattice is essentially to fix the values of the 
coefficients {^m}. In the following, wc attempt to expand the solution to the spin distortions of an x-directed AFM 
bond near the origin in the basis {cos(m6')/r™}. 

To start, we expect that the coefficient Ai must dominate the others since cos{9)/r is the mode of lowest energy. 
Further, convergence at the bond sites requires that {Am) faster than 2~™ as to ^ oo. Thus, it is meaningful 
to treat the expansion 0("^(r, 6), consisting of the first n terms of Eq. ( p^ as an approximate solution. We can then 
apply the local equilibrium condition at n sites around the bond to determine the n coefficients. 

For concreteness, consider the four term expansion 

, , ,(A\, „^ . cos 9 , cos 30 , cos5(? , cos7(? 

0(r, 9) = 0W 9)=Ai + A^^^ + A,^^ + . (12) 

To solve for its four coefficients we require the 9x4 transformation matrix 
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(13) 



(generated using a symbolic algebra computer package) and the matrix 
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(14) 



of linearized equilibrium conditions at sites 1 through 4. Defining the row vector a''" = [j4i, ^43, Ay] the deter- 
mination of {Am} is equivalent to solving the homogeneous system of equations MTa. = 0. The requirement that 
det(AfT) = yields a critical value Ac = 1.0113 (very close to the true value Ac = 1) for which 



a = [0.5987, 0.0742, 0.1960, -0.0552] x 
That is to say, the best four term expansion reads 

cos 39 



A, 



A. 



COS 59 



A. 



COS 79 



(15) 



(16) 



m— 3,5,7 



with 



Ai 



-0.5987(/)i, A3 = +0.074201,^5 = +0.1960(/)i, A7 = -0.0552(^i 
The coefficients of 0'-"-' (r, 9) for n being increased from 1 to 4 are presented in Table |. 

TABLE I. Spin distortion amplitudes of (r, 6) for n = 1, 2, 3, 4. 
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What this calculation provides that the other does not is the value of the multiplicative factor dAi/d(l>i = Ai/(f)i ^ 
0.6 relating the magnitude of the spin distortions at the bond sites to the magnitude of the dipole moment associated 
with the bond itself. 

We have shown that the spin distortions are given everywhere by 

,, ^. , cos^ , cos 30 , cos 50 

cPir, e)=A,^+A,—^+A,—^ + ... (18) 

As we have seen, however, this expression is unwieldy in that it requires the application of infinitely many local 
equilibrium conditions to fully determine the coefficients {^to}- Even to calculate the coefficients of a finite series 
expansion of several terms is computationally expensive. Ideally, what wo would like to have is a solution dependent 
on a single parameter whose value is determined by applying a single boundary condition at the bond itself. Here 
we now outline such a method. We find that a simple assumption on the distribution of modes can provide this very 
result. 

We proceed by assuming that the spectrum of modes can be modeled by 

2/1, 

or, in somewhat simplified notation, 

= (±) fc^ (20) 

where the index k is taken to be odd and the sign is taken alternately positive and negative. This form falls off just 
fast enough to make the series converge — besides this seemingly naive reason, we appeal to its success (described in 
detail below) to justify its usage. 
Under this assumption 

,, , COS0 , cos 30 , cos 50 , 

<^(r, 0)=A,--+ A,—^ + + . . . (21) 

, / cos 1 cos 30 1 cos 50 
= 2Ai — :T77rTT + 



2r 3 (2r)3 5 (2r)5 

Now, the magnitude of the dipole moment in terms of the distortion at the bond site follows immediately from solving 
(/)(l/2, 0) = (^1 self-consistently. We find that Ai = and hence 

,/ , 4 v-^, , , cosfc0 

^(r,0)=^,-B±)^ (22) 

since 

-0(1/2,0) = 0ii(^l-i + i-...) (23) 

which is identically equal to The spin distortions at the remaining sites in the immediate neighbourhood of the 
bond are as follows: 

02 = ^(V5/2, arctan(2)) = 0.295167235301 , 04 = (/'(3/2, 0) = 0.409665529401 (24) 

(In fact, one may prove the identity 202 + 04 = 0i, which we shall use later on in this paper.) 
Notice that as r becomes large, we get 

4 / 1 \ 4 cos 

0(r, 0) — > 01 — cos arctan ( — ] 6i (25) 

TT V 2r / TT 2r 

, 2cos0 



= (pi 

TT r 

so that the solution retains its familiar long range behaviour. That is 



•/-(r) = ^ (26) 
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but now with 



p=\p\ = -n 

TT 



(27) 



What remains is to determine the parameter As promised, the bond furnishes a single boundary condition in 
the form of the equihbrium condition apphed at either of the bond sites: 



Xsm{2(f>i) + 2sin((/)i — </'2) + sin( 



i) = 0, 



(28) 



This is an imphcit equation for 0i, and thus for all of the spin distortions as a function of A. Its solution is plotted 
in Fig. 0. 
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FIG. 2. The distortion magnitude as a function of the frustrating AFM bond strength A. The solid line represents the 
prediction of our analytical theory, and the broken line is the prediction based on the solution of Kovalev and Bogdan [ p3[ . 
Open squares are numerically generated data points from large lattices (from 40 x 39 up to 80 x 79 lattices) with converged 
numerical solutions for this quantity. 



Moreover, the linearized equation gives 

- Ac20i + 2(</.i - 02) + (01 - 04) = (29) 



which can be solved explicitly: 



3_202+04^^_ (30) 
2 201 ^ ' 



That is, our ansatz correctly reproduces the exact critical value of Ac ! 

A comparison of the solution Eq. ( ^2|) to numerical simulations is presented in Fig. |3[ clearly, the agreement is 
excellent, providing the most direct support for our ansatz. 
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FIG. 3. The local spin distortions as a function of distance along the a;-axis for an a:-directed AFM bond. The solid line 
represents the function 0(r, 0) described in the text. Open squares are converged numerically generated data points from 
simulations performed on the same size lattices as in Fig. 2. The broken line is the (^'^•'(r, 6^) solution of Kovalev and Bogdan 
[p3[ which fails to predict quantitatively the correct dipole moment of the long-range distortions. 



C. Energy functional 

We now calculate the total energy stored in the spin distortions induced by a single AFM bond. That these 
distortions arc both small and slowly varying in position away from the core allows us to convert the sum of the 
energy contributions into an integral of the energy density (V(/))^. A complete derivation is provided in Appendix A, 
and we summarize the results below. 

To begin, the Hamiltonian is approximated to second order everywhere except across the AFM bond itself: 

HK-Y,Jr] + \JY^{<I^^ - (1^3? - 2AJsin2 01 (31) 

Of course, the term — ^ represents the total energy of the system in the absence of spin distortions so that the 
energy from the distortions alone is given by the latter two terms. They, in turn, can be expanded using the continuum 
approximation 

Eiist ^\jY. - ^^+-^^ + l^i'^i^^ - '^4)' + 4((/.i - 02)') - 2AJsin2 0i 




where M is the xy-plane excluding a small region about the bond centre. 

An explicit evaluation of the energy using the solution of the previous subsection is presented in Appendix A, 
wherein the full effect of the core region is accounted for. We find 

Sdist = 2j(02_Asin2 0i) . (33) 

This result is compared to numerical simulations in Fig. ^, and again, excellent agreement between our numerical 
solutions and our analytical work is found. 

So, now we carry on to the examination of the many-bond problem, having an excellent solution to both the core 
and far-field distortion patterns of the single-bond problem, as well as an accurate energy functional for an isolated 
frustrating bond. 
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FIG. 4. Energy of a single bond (in units of J) as a function of the relative AFM bond strength A. The solid line represents 
the prediction of theory, and open squares are converged numerically generated data points from simulations performed on the 
same size lattices as in Fig. 2. 



After dealing with a single bond in isolation, the next step toward treating a non-zero density of bonds is to 
determine how bonds interact with one another. In this subsection we consider the problem of two bonds. 

Suppose that there is an AFM bond, call it A, between the and 1 sites. Then suppose that another similarly 
directed bond, call it B, is placed between the s and s + 1 sites and that a vector R making an angle <i> with the 
X-axis connects the two bond centres, as in Fig. 0. 



III. INTERACTING FRUSTRATING BONDS: 







1 



s 




FIG. 5. The geometry and labelling used for two parallel bonds separated by a vector R. 



We expect that the energy can be parametrized by two variables 



a = 



4>i - 4>o 

2 



0S+1 - 4>s 

2 



(34) 



and that it can be written in the form 



E2{a, (3) = E{a) + E{I3) + Jga(3 




(35) 
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where g = (;(R) is a real- valued function of the separation and relative orientation of the bonds and J is a small 
correction representing the 4**^ order contribution to energy from the spin distortions which were neglected in Eq. (|3^). 

The requirement that the spin texture remain invariant up to a global sign change under interchange of a and (3 
reduces the two parameter energy expression (|35| ) to one of two one parameter expressions or E^ corresponding 
to the symmetric [a = j3) and the antisymmetric (a = —(3) state. 

Case 1: (j) — a ~ —(3. In this case, the dipoles associated with the bond are anti-aligned. The total energy is 

E:,{cp)^E2{^,-<f)=Aj[{\-\--^g)<l,^ + \{l-5]x4>'^^ < 2E{4>). (36) 
Minimization with respect to (j) gives 



2 V3 



^=^|{\ + \9-l)|(^^-5^\. (37) 

This implies that the critical value of A at which the canted ground state first appears is lower than it is for a single 
AFM bond Ac = 1. 

Case 2: (p = a = f3. This represents a higher energy metastable state characterized by aligned dipoles. The energy 
for this configuration is 

i?2+(0)=i?2(0,0)-4j(^(l-A+i5)02 + i(^^_<5^A(^4^ > 2E{<P) (38) 
with distortion magnitude 



= W(A--.9-l)/ --<5 A. (39) 



Now the critical value of the coupling constant required for a distorted ground state exceeds Ac = 1. 

Case 1 is of particular interest since it yields the ground state energy E2 of the two bond system. Further, by 
re-expressing that energy in terms of the energy of a single bond (extracted from Eq. (|3^ ) in the limit A — s- Ac -I- 0"'") 
we can determine the energy of interaction between the two bonds. 

In general, since we expect g to be small, we can write 
However, as A ^ 1, we obtain 

a weak, long-range interaction with a higher power law, a consequence of the fact that the dipolar distortions do not 
pre-exist in the unperturbed medium at A = 1. Such an interaction is analogous to the Van der Waals interaction 
between thermally fluctuating dipoles. 

The function g expresses the functional dependence of the interaction energy on the geometrical configuration of 
the bonds. We have yet to determine its exact form. All we can say now is that 

hm g{R) = (43) 

it— »oo 

which simply formalizes our expectation that two bonds must be non-interacting at infinite separation. 

The long range spin distortions arising from each of the bonds A and B with dipole moments p"^ and is given 

by 

A B 

</.-4(r) = P-^ and </,S(r) = . (44) 
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The linearity of the Laplacian imphes that the field equations admit a superposition principle. Therefore, we take the 
total spin distortion at each point to be 



(45) 



where and (p^ are the spin distortions from each bond in the absence of the other. 
The energy in the distortions away from the cores goes as 



/ (V</>fdV=/ (V(</>^ + 0^)2)dV (46) 

Jm2 Jm2 



Jm2 Jm2 



I M2 J M2 

JM2 

where M2 (= R\L)(:(0)\-De(R)) is the xy-plane with disks removed about the bond centres. The last term in this 
expression vanishes identically, which indicates that the core of each bond interacts with the long-range spin distortions 
of the other and justifies a rather involved calculation of the interaction energy which we have relegated to Appendix B. 

What we find is that E'int has the form of a magnetic dipole interaction. Further, although in the preceding 
discussion we considered only parallel bonds, it is simple to show that these results hold more generally. Thus, for 
two bonds which are parallel or perpendicular we have 

Sint = {2(p^ • R)(p^ • R) - • } . (47) 

Finally, we can work backwards to find ^(R). For parallel bonds, 

A(2/3 - 5) -kB? 

That is 



B., = J% {2(p- . RXp" . R) - . p"} = ±,^^-1,0032* . (48) 



9(R) = -Hj':os2*. (49) 
We note that the identical calculation for two perpendicular bonds gives 

5(R) = ^sin2$. (50) 




FIG. 6. The ground state energy of two parallel bonds, for A = 1.1, at various separations along the x-axis, where the solid 
line is our prediction and the open squares are converged numerically generated data points from simulations performed on the 
same size lattices as in Fig. 2. 
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FIG. 7. Energy of two parallel bonds, A = 2, at various separations along the a;-axis, where the solid lines and open squares 
are converged numerically generated data points from simulation performed on the same size lattices as in Fig. 2. The upper 
curve corresponds to the metastable state with dipoles similarly directed (-^^), while the lower energy curve corresponds to 
oppositely directed dipoles (—><—)• 



IV. CLUSTERS IN THE CLUSTER SPIN GLASS PHASE: 



We should ask whether the results we have obtained so far for the one and two bond problems can be generalized 
to allow us to tackle the problem of a lattice frustrated by the presence of any number of arbitrarily placed AFM 
bonds. It should be clear that, in general, even for relatively few bonds, the induced spin distortions will be very 
complicated and the energy surface characterized by many closely spaced, low-lying states. In such a case one must 
resort to sophisticated computer algorithms to numerically generate the ground state spin texture, and the qualitative 
analysis of the cluster spin-glass phase from such work has been analyzed elsewhere p5|,p[ . 

In contrast to that work, here we note that there are configurations of suitably high symmetry for which we can 
confidently treat the spins as planar and even solve analytically for the spin distortions and energies of all the possible 
states. The smallest such configuration is the square cluster of parallel bonds. 

By a cluster we imply a collection of bonds arranged in some local region on the lattice. That collection, call it 
C, can be thought of as a set of dipole-position pairs: {(p", r")}^^^ • Given a high symmetry cluster, for which a 
spin-planar ground state is justified, the spin distortions away from the cores of the bonds are given by 

aec ' ' 

which is the solution (unique for the required symmetry) to the equation 

V20(r) = -27rV • ^ p"(5(r - r") . (52) 

aec 

The total interaction energy of the cluster can be written as the sum of all pairwise interactions 
p \- . 277 ,^ p"-(r"-r^)p^.(r"-r^) 

Thus, for instance, the L x L square cluster of parallel bonds given by 

ri = (0,0), r2 = (L,0), r3 = (0,L), r'* = (i,L), p" ^ ±px for a e C = {1, 2, 3, 4}, (54) 
has a rather simple interaction energy. There are four possibilities depending on the the orientation of each dipole: 
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Eint = -J— p\ 0,0, +J—p^ 



(55) 



In practice, however, the degeneracy of the middle two states is hfted by higher order terms in the interaction energy. 
Indeed, the sphtting observed in numerical simulations enables us to list the four distinguished states in order of 
decreasing energy. 



first excited state 



ground state 



(56) 



The first excited state consists of four similarly directed dipoles. Destructive interference inside the cluster gives 
zero net distortion, but outside the cluster the dipoles add constructively so that the cluster acts like a single unit 
with a much stronger moment. This results in very strong long range spin distortions. Far enough from the bond, 
the spin distortions are given by 



(r) 



4p • r 



(57) 



Since the long range spin distortions mediate the interaction between bonds, we expect a cluster of this kind to 
strongly couple to other bonds in the lattice. 

In the ground state, for which the spin distortion pattern is shown in Fig. ^ we have the opposite case: internally, 
the dipoles add constructively to give large distortions whereas outside they cancel to give very small ones. The 
absence of long range spin distortions implies that these clusters can only weakly interact with other bonds. Most 
interesting, though, is that the internal spins are uniformly oriented but differently ordered from those spins outside 
the cluster. 
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FIG. 8. The ground state spin texture of the 4x4 square cluster for four frustrating bonds. The pink spins are bond sites 
at which the frustrating bonds are placed. For clarity, the spins within the locally ordered cluster are coloured red. 



Numerical solutions of the energies of this square cluster are shown in Fig. ^ The solid curves are our analytical 
results (apart from a constant times the single-bond energy (that is straightforward to calculate)), and provides strong 
support for their usage. In the ground state this figure makes clear that the binding energy of the cluster can be quite 
large, especially for small cluster sizes. Further, square clusters tend to settle into states with strong internal binding 
and which interact only weakly with other bonds. That is to say, a square cluster is a locally ordered domain whose 
local order parameter is non-coUinear with spins in the rest of the lattice. However, this behaviour is a strong 
bmction of cluster size. This construction, and the energy plot of Fig. ^, demonstrates that such clusterings of spins 
(in regions that are not too large) are stable. 
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FIG. 9. Energy of the ground (open squares and lower curve) and first excited (open diamonds and upper curve) states of the 
four frustrating bond arrangement described in the text. The open squares and diamonds are converged numerically generated 
data points from simulations performed on the same size lattices as in Fig. 2. 



The relation of such clusters to stripes, or in the case of spin modulations, to the physics associated with the 
appearance of domain walls, can be demonstrated by considering the interface between such clusters, and to this end 
we have analyzed a six-bond cluster given by 

ri = (0,0), r2 = (L,0), r3 = (2L,0), v^ = {0,L), r' = {L,L), r^^{2L,L), (58) 
p" = ±pi for a e C = {1,2,3,4,5,6} . 

Following the above analysis for four bonds, for the six-bond situation there are 64 possible choices of the dipole 
moments' orientations, and these states have 15 different energies. The lowest energy configuration corresponds to 




(59) 



and has a dipole-pair interaction energy of (— 1037rJ/10) (p/L)'^ (which is noticeably lower than the first excited state, 
which has an energy of (— 387rJ/10) {p/L)^). 

The ground state spin texture for this location of the six frustrating bonds is shown in Fig. From this figure one 
can see the important result that this arrangement of spins is exactly what would expect if each 4x4 cluster within 
the 6-bond cluster was in its respective ground state. Thus, between these 4x4 clusters one obtains a domain wall, 
of width one lattice spacing, over which the local magnetic order parameter is rotated. Numerical evidence suggestive 
of this type of domain wall was discussed at length in Ref. for the case of a random distribution (and orientation) 
of a non-zero density of frustrating spin interactions, and it is clear that the same physics is at work in these two 
situations: spin twists. 
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FIG. 10. Ground state of the six frustrating bonds showing the domain wall between two clusters of locally ordered spins. 
As in Fig. 8, the spins of the bond sites are coloured pink, while the two clusters are coloured red and green, respectively. 



V. CONCLUSIONS: 

The above formalism has provided a detailed analytical theory of the spin distortions generated by frustrating 
bonds, and of the interactions mediated by the spin background between them. We have critiqued its validity by 
comparing to numerical solutions, and have found excellent agreement. Then, by focusing on highly symmetrical 
distributions of frustrating bonds, reminiscent of local regions of bonds in the multiply doped state, we have used 
this theory to verify the existence of locally ordered magnetic clusters. Most importantly, our work shows that these 
clusters are stable. These are the clusters envisioned to exist in the so-called cluster spin-glass phase of LSCO 
and Yi_,CaxBa2Cu3)6 @. 

It is to be stressed that it is believed that mobile holes produce the same kinds of spin distortions that are 
produced by the frustrating bonds discussed in this paper ]l^ , p6| . Thus, we believe that our results support previous 
suggestions p"^-p^ that spin twists and distortions are part of the competing interactions that might lead to rivers 
of charge appearing as low-energy fluctuations in the doped cuprate systems. 
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APPENDIX A: SINGLE BOND ENERGY 



The question of the energy stored in the spin distortions can be answered by expanding the Hamihonian as foUows: 

H = — cos (0i — + A J COS (01 — ^o) 

«-J^' +AJcos(20i) 

+ ^'^Il'(<^'' - -^j)" - 2AJsin2 01 (Al) 

(ij) (ij) 

Of course, the term — ^ J^j represents the energy intrinsic to the lattice. Therefore, the energy from the spin 
distortions alone is given by 

i J^'(0, - 0j )2 - 2AJsin2 01 . (A2) 
This in turn can be expanded using the continuum approximation 



Edist ~J{ \ I (V0)"d2r + (01 - 04)2 + 2(01 - 02)2 - 2Asin2 0i 

J M 



0V0 • ds + (01 -04)2 + 2(01 -02)2- 2Asin2 01 (A3) 



where M — TV\D^{0) is the xy-plane excluding a small region about the bond centre. The integral term of this 
expression must be evaluated. Since the core of the bond occupies a unit disk at the origin, the appropriate value for 
e is 1. Thus, one finds 

i^dist ^j[\Y.^^l + ("^1 - '^4)2 + 2(01 - 02)2 - 2Asin2 0i| . (A4) 

^ k ^ 

Using the results for the mode characterization of Eq. (|l^) gives 

k k ^ ^ ^ 

and the identity 202 + 04 — (/>i allows us to write 

(01 - 04)2 + 2(01 - 02)2 = 02 + 202 + 02 . (A6) 

Thus the energy in the spin distortions is 

Eiist - J ( f - In ( I ) + 1 ) 02 + 202 + 02 - 2A sin2 0i I (A7) 



which can be evaluated using Eq. ( p^ to give 

£;dist -2J(02- Asin2 0i) . (A8) 
We stress that this expression includes the energy from the core of the spin distortion field. 
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APPENDIX B: INTERACTION ENERGY 



In accordance with Fig. ||, we write the spin distortion at the site from the bond B as 



p^-(R+^x) 



p^.(R 



R X 1 



i?2 4i?2 



(Bl) 



A similar expression can be written down for the other bond site. Then, by superposition, the difference between the 
net distortions at sites and 1 is 



1 



i)^ic^^-cbt)-—ip^--R){p^-x) + —p^-x + — 



i?4 



1 



i?2 



(B2) 



Then we make use of the relationship between the spin distortions at the bond sites and the magnitude of the dipole 
moment: 



Therefore 



so that 



(</>o - </>i)^ = (00^ + c^tr -^Ap^- R)(p^ • R) + Sp^ • 



1 



1 



< — ^1 



2" ' 2 \R^' 

and thus, by symmetry, the total energy is 



An 
R^ 



(p-4.R)(p^-R) 



E = Eip^) + Eip^) + J{^p^ . - g(p^ . R)(p5 . R)} 



We conclude that the interaction energy is given by 



.2tt 



i?int(p",p") = |-^(p^ . R)(p^ . R) - p^ . p^ 

(Repeating the above calculation for two perpendicular bonds produces the identical result.) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 
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